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SUMMARY  OF  RESULTS 


Several  algorithms  for  self  cohering  an  airborne  distributed  array 
on  radar  clutter  have  been  tested  and  their  performances  have  been 
statistically  evaluated.  Also,  several  new  algorithms  have  been 
developed  throughout  the  course  of  the  study.  The  new  algorithms 
offer, in  general,  improved  performance  and  wider  applicability. 

We  obtained  airborne  clutter  data  from  two  sources  (Task  £)  :  the 
Naval  Research  Laboratory  (NRL)  and  the  Environmental  Research 
Insitute  of  Michigan  (ERIM).  We  identified  and  obtained  specific 
runs  of  the  8-Element  NRL  airborne  array.  These  covered  land  as 
well  as  sea  clutter,  different  radar  look  angles,  and  different 
pulse  repetition  rates.  We  also  obtained  several  runs  of  ERIM  SAR 
data.  We  adapted  the  use  of  SAR  data  to  permit  testing  of  the 
"Minimum  Variance"  and  “Robust"  algorithms. 

We  tested  and  statistically  evaluated  the  performance  of  the 
Minimum  Variance  algorithm  (Task  1)  using  the  NRL  and  ERIM 
clutter  data.  The  algorithm  was  successful  against  land  clutter 
where  it  was  used  successfully  for  finding  targets  of  opportunity 
that  allowed  retrod irect ive  beamforming.  The  algorithm  was 
unreliable  (as  expected)  against  sea  clutter. 

We  also  tested  Doppler  Beamsharpening  techniques  (Task.  3)  against 
the  NRL  land  clutter  _  data.  Doppler  Beamsharpening  proved  to 
be  highly  effective  in  increasing  the  chance  of  finding  targets  of 
opportunity  that  allow  high  quality  retrod irect ive  beamforming. 

We  developed  several  new  algorithms  (Task  4)  during  this  study. 
We  developed  and  tested  the  Robust  Minimum  Variance  Algorithm  to 
provide  reliable  beamforming  on  a  collection  of  inferior 
bearnformers  (targets  of  opportunity).  We  also  developed  and 
tested  the  Spatial  Correlation  Algorithm  which  is  applicable  to 
sea  as  well  as  to  land  clutter.  It  allows  also  for  local 
processing  within  the  aperture,  in  contrast  to  transmitting 
information  from  each  module  to  a  centralized  rigid  processor, 
which  makes  it  attractive  for  a  cost  effective  large  aperture 
space  based  radar.  Finally,  we  developed  a  technique  of  Subarray 
Processing  for  applying  the  above  new  algorithms  on  very  large 
airborne  or  spaceborne  radar  systems. 


INTRODUCTION 


The  objective  of  this  effort  is  to  investigate  techniques  that  can 
be  used  to  mount  airborne  phased  array  radar  antennas  on  arbitrary 
portions  of  the  airframe.  Current  airborne  phased  arrays  cannot 
be  so  mounted  because  the  nonrigidity  of  the  surface  would 
severely  degrade  the  phase  coherence  necessary  for  effective 
phased  array  performance.  The  ability  to  distribute  the  elements 
of  a  phased  array  over  the  entire  aircraft  will  greatly  enhance 
the  ability  of  airborne  radar  systems  to  function  in  sophisticated 
battle  environments. 

The  performance  of  an  airborne  radar  antenna  array  is  limited  in 
several  highly  important  respects  by  the  fact  that  the  array  can 
be  mounted  only  on  certain  sections  of  the  aircraft.  This 
localisation  of  the  antenna  results  in  incomplete  surveillance 
coverage  and  restricted  freedom  of  communication  with  all  desired 
stations.  The  restricted  antenna  aperture  also  limits  the  usable 
detection  range  of  the  radar  and  the  angular  resolution  of  the 
system,  and  makes  the  radar  more  susceptible  to  mainbeam  jamming. 
These  limitations  could  be  greatly  reduced  if  arrays  could  be 
installed  in  any  arbitrary  portion  of  the  airframe.  Air  Force 
supported  work  in  the  past  five  years  has  shown  that  there  is 
promise  in  designing  radar  antenna  arrays  to  make  use  of  arbitrary 
portions  of  the  airframe.  The  key  idea  behind  the  use  of  such 
arrays  is  to  make  use  of  a  portion  of  the  scattered  return  of  the 
radar  signal  to  self-cohere  the  phases  of  the  array  elements  and 
so  compensate  for  the  distortion  in  phase  caused  by  the 
nonrigidity  of  the  skin  and  airframe  that  would  otherwise  severly 
degrade  the  antenna  performance.  Not  ail  scattered  radar  returns 
are  useful  for  this  purpose,  however,  and  it  is  important  to 
investigate  the  statistics  of  realistic  radar  clutter  returns  to 
see  if  systems  can  be  designed  to  utilise  these  returns  for 
sel f-synchroni z ing.  The  objective  of  this  work  is  to  identify 
existing  sources  of  useful  data  and,  using  these  data,  to  test  and 
compare  various  algorithms  for  synchronising  the  array. 


CLASSIFICATIONS  OF  SELF-COHERING  ALGORITHMS 


Table  3.1  gives  a  classi f icat ion  of  the  self-cohering  algorithms 
under  consideration  (including  the  newly  developed  algorithms). 
The  algorithms  are  classified  both  according  to  the  type  of 
processing  and  according  to  the  environment. 

Regarding  the  type  of  processing,  we  classify  the  algorithms 
according  to  whether  or  not  Doppler  processing  is  used.  It  is 
worthwhile,  however,  to  note  that  Doppler  and  non-Doppler 
techniques  can  be  combined  in  some  cases  to  come  up  with  superior 
algorithms.  For  example,  and  as  will  become  clear  later  on,  the 
Doppler  Bearnsharpening  techniques  can  be  si  gni  f  icant  ly  enhanced  by 
coupling  the  Doppler  time— domain  processing  with  the  Robust  space 
(or  range  bin) -domain  processing. 

With  respect  to  the  environment ,  we  distinguish  between  the 
situation  where  the  reradiating  scatterers  are  rigid,  which  is 
found  only  in  land  clutter,  and  the  situation  where  the  scatterers 
undergo  appreciable  movement  (in  units  of  wavelength)  due  to  wind 
as  found  in  vegetation  and  sea  clutter. 

TABLE  3. 1  CLASSIFICATION  OF  SELF-COHERING  ALGORITHMS 
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4.  SIUDY^RESULTS 

In  this  section  we  give  a  detailed  report  of  study  results.  The 
section  is  organized  according  to  contract  tasks.  For  the  sake  of 
continuity,  lengthy  details  and  long  mathematical  deviations  are 
defered  to  appendices. 


4.1  TASK  1:  MINIMUM  VARIANCE  ALGORITHM 

Figure  4.1  shows  a  high  angular  resolution  bistatic  radar 
system  that  consists  of  a  pulsed  transmitter  and  a  large 
distributed  thinned  array  of  receivers.  In  principal,  the 
system  can  be  ground  based  or  airborne.  For  a  ground  based 
system  the  size  of  the  receiving  array  can  be  increased  to  the 
extent  that  angular  resolving  powers  comparable  to  common 
optical  instruments  such  as  cameras  and  small  telescopes  can 
be  achieved.  This  system  has  been  introduced  in  the 
literature  111,  C21  as  a  "radio  camera"  capable  of  "taking 
pictures"of  microwave  scattering  objects.  For  an  airborne 
application,  the  size  of  the  receiving  array  will  be  limited 
only  by  the  size  of  the  airframe.  Yet,  and  as  mentioned 
earlier,  an  airborne  phased  array  radar  flush  mounted  on 
arbitrary  portions  of  the  airframe  offers  many  advantages  over 
existing  phased  array  systems  which  tend  to  be  much  smaller  in 
size  in  order  to  maintain  a  rigid  structure. 

Such  a  large  airborne  antenna  system  has  to  be  self-adaptive, 
i.e. ,  the  beam  must  be  formed  by  a  self-cohering  process. 
Otherwise,  the  nonrigidity  of  the  surface  would  severly 
degrade  the  phase  coherence  necessary  for  effective  phased 
array  performance.  The  basic  idea  is  to  find  a  suitable 
reference  source,  use  adaptive  ret  rod  i  rect  i  ve  bearnf  orrn  l  ng 
techniques  to  focus  the  large  distorted  array  on  it,  and  then 
move  the  narrow  focused  beam  in  range  arid  angle  to  a  target 
area  and  scan  the  beam  across  the  target. 

It  is  shown  in  C13,  C23  that  a  suitable  reference  should 
approximate  a  point  source,  strong  enough  to  overcome  clutter 
(about  6  dB  higher).  That  is,  it  should  have  the  right 
combination  of  small  physical  size  (in  terms  of  wavelength) 
and  strong  echo.  In  a  ground  based  system,  the  phase 
synchroni z ing  source  can  be  a  corner  reflector,  a  beacon,  or  a 
target  of  opportunity  (TOO).  In  an  airborne  system,  the  basic 
algorithm  should  rely  only  on  finding  suitable  TOO’ s.  The 
minimum  variance  algorithm  was  developed  to  do  exactly  that 
C23. 

4.1.1  PROCEDURE 
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Figure  4.1  describes  the  procedure  explicitly.  A 
ground-to-ground  imaging  situation  is  portrayed.  On  the  left 
are  a  radar  transmitter  and  receiver.  The  target  region  to  be 
imaged  is  on  the  right.  A  transmitted  pulse  illuminates  the 
target  area  and  a  passive  phase-synchronizing  target,  which  is 


sketched  as  a  corner  reflector.  The  corner  reflector  is  shown 
closer  to  the  radar  than  the  target  region.  Each  antenna 
element  receives  an  echo  trace  that  contains  clutter,  the 
beamforming  corner  reflector,  and  the  target.  The  echo 
waveform  from  each  element  is  coherently  demodulated,  sampled 
in  range,  and  stored  in  the  format  shown  for  image  formation. 
To  produce  the  image  the  processor  must  perform  several 
functions:  it  must  find  the  range  element  that  contains  the 
phase  synchronizing  source,  self-focus  a  beam  on  it,  refocus 
the  beam  to  the  target  range,  and  scan  the  beam  in  azimuth 
across  the  target.  When  these  operations  are  completed,  the 
output  of  the  system  is  a  one-dimensional  scan  or  image  of  the 
target.  A  two-dimensional  (range,  angle)  image  is  built  up 
from  a  succession  of  angle  scans. 


THE  RADIO  CAMERA 
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FIGURE  4.1  RADIO  CAMERA  DATA  FORMAT  AND  PROCEDURES  (FROM  CEO  ) 


The  range  trace  delivered  by  each  antenna  element  is  stored  as 
a  row  of  complex  numbers.  Successive  rows  correspond  to 
successive  elements  in  the  array.  The  position  of  a  sample  in 
a  row  is  proportional  to  range  and  designates  the  range  bin. 
The  position  of  a  sample  in  a  column  designates  the  array 
element  number  and  is  monotonic  with  but  not  necessarily 
proportional  to  element  position  because  the  array  may  be 
distorted. 
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sampling  and  storage 
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range  from  which  the  phase  synchronizing  reference  signal  is 
obtained.  The  radiation  field  from  a  point  source  in  free 
space  would  be  nearly  constant  in  amplitude  across  the  array, 
whereas  the  phase  differences  would  disclose  the  differential 
distances  (modulo  wavelength)  to  the  source  as  well  as  the 
differential  phase  shifts  through  the  receiving  elements.  A 
simple  test  to  find  the  target  with  a  radiation  field  that 
most  closely  approximates  that  of  a  point  source  is  to  measure 
the  normalized  echo  amplitude  variance  at  each  range.  RQ  is 
that  range  for  which  the  test  value  is  minimum. 


Let  us  denote  the  echo  amplitude  frorn__  the  ith  range  bin 
received  by  the  nth  element  ftjn  .  Its  mean  ftj  =  <1/N )  <ancl 

mean  square  ft*  =  (1/N)  £ ,  A;*  are  calculated,  where  „  N  is  t_he^ 
number  jOf  antenna  elements.  Its  variance  is  (1/N)  (Ain-  ft,-  ) 

=  ft{  -  ft,.  .  The  variance  is  normalized  to  ft*-  to  remove 

the  effect  of  target  strength  and  range, from  the  calculations. 
The  normalized  variance  then  becomes  &;■  -  1  -  ft.  /ft1;  .  The  rule 
for  finding  R„  ,  then, is  to  calculate  07  for  all  range  cells 
and  search  for  the  smallest  value.  The  signal  sequence  from 
that  range  bin  becomes  the  reference  signal.  ftppendix  ft  gives 
a  detailed  description  of  the  adaptive  beamforming  and 
scanning  operations  that  result  in  the  high  resolution 
microwave  image. 


It  is  worthwhile  to  list  the  following  general  features  of  the 
algorithm: 


Simple  and  fast. 

Beamforming  angle  can  be  anywhere  in  the  illuminated 
sector  (especially  for  an  airborne  application  where  we 
entirely  depend  on  finding  TOO’ 5). 

Requires  individual  amplitudes  from  every  element  to  be 
communicated  to  a  central  processor. 
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4.1.£  BEftMFORMER  QUALITY 


The  loss  in  mainlobe  gain,  in  dB,  due  to  all  the  random  phase 
errors  across  a  linear  array  (having  a  large  number  of 
elements)  is  shown  in  Cl  3,  C23  to  be 


AG  -a  4.3 


(4.  1 ) 


where  is  the  variance  of  the  phase  errors  i assumed  to  be 
normally  distributed)  in  square  radians.  The  relationship  1.3 
plotted  in  Figure  4.  3.  The  loss  is  minor  <  <1  dB)  for  an  RMS 
phase  error  of  about  1/3  rad.  The  loss  grows  rapidly  as  the 
phase  error  grows.  A  1  rad  RMS  phase  error  resuits  in  ■‘t.  3  aB 
loss. 


FIGURE  4.3  DECIBEL  LOSS  IN  GAIN  DUE  TO  GAUSSIAN  ERRORS  (FROM  L21 ) 

In  our  statistical  evaluation  of  the  Minimum  Variance 
algorithm  the  standard  deviation  of  the  residual  phase 
errors  across  the  array  at  the  ueamformmg  angie  _  3  used  3 
our  measure  of  quality  of  candidate  beamforming  sources.  We 
were  able  to  estimate  because  the  radar  clutter  data  that 
we  examined  were  collected  by  well  designed,  undistorted 
arrays  with  high  clutter  to  noise  ratios  (the  NRL  3-element 
array  and  ERIM  3AR).  This  means  that  the  measured  phase 
across  the  array  (for  a  particular  range  ceil)  is  almost 
entirely  due  to  the  illuminated  clutter  patch.  In  a  real- life 
system  of  a  distributed  airborne  array,  array  distortion  will 
obscure  the  received  phasefront  and  precludes  the  estimation 
of  . 

A  good  beamforrner  m  the  far'  field  of  a  linear'  array  should 
produce  an  almost  linear  phase  across  it  with  the  possibility 
of  discontinuities  at  iTT  due  to  the  modulo  E  nature  of  pnase 
measurement s.  1  he  slope  of  the  underlying  straight  line 
discloses  the  angular  position  of  the  dominant  scatterer 
within  the  illuminated  patch.  Retrod  meet ive  beamforming 
forms  a  beam  focused  on  that  scatterer'  leaving  the  difference 


between  the  underlying  straight  line  and  the  measured  phase  as 
a  residual  phase  error  pattern  across  the  array.  To  estimate 
its  standard  deviation,  we  first  unwrap  the  phase  (if 
necessary)  by  adding  or  subtracting  the  appropriate  multiple 
of  £  TT  (that  removes  any  discontinuities  in  the  measured  phase). 
A  linear  regression  then  follows  and  the  difference 

S  (X)  between  the  unwrapped  phase  and  the  regression  line  is 
obtained.  Figure  4.3  illustrates  these  steps.  An  estimate  of 
is  obtained  as 

A  a/  1/2 

=  (  h  l  )  (4-£) 

H9* 

where  N  is  the  total  number  of  antenna  array  elements. 

,  as  given  by  <4.£),  together  with  (4.1),  will  be  our 
quality  measure. 


FIGURE  4.3  PHASE  UNWRAPPING  AND  LINEAR  REGRESSION 


§I0IISIICS_OF  crz<p 

In  both  the  NRL  Q-elernent  array  data  and  ERIM  SAR  data  we 
applied  the  Minimum  Variance  Algorithm  repeatedly  on  different 
runs  examining  each  time  the  range  bin  that  has  the  minimum 
amplitude  variance  and  estimated  as  described  in  the 
previous  section.  We  also  examined  other  range  bins  that  had 
low  amplitude  variance.  Although  we  consistently  found  that 
the  range  bin  of  the  minimum  amplitude  variance  had  a  low  or  , 


-9— 


it  was  net  necessarily  the  best  range  bin.  Other  range  bins 
(of  low  amplitude  variance)  were  often  better  in  the  sense 
that  they  gave  lower  •  We  concluded  that  it  is  more 

meaningful  to  examine  the  statistics  of  for  these  range 

bins  that  produced  amplitude  distributions  of  low 
variance.  The  procedure  consists  of  setting  a  threshold  on 
the  normalized  standard  deviation  of  the  amplitude,  ^  , 

and  examining  the  phases  produced  by  all  the  range  bins  that 
are  below  this  amplitude  variation  threshold.  We  then  sort 
the  resulting  set  <  a#.  \  rA<  <  ^  >  in  an  ascending  order  and 
produce  an  estimate  of  the  cumulative  distribution  function 
(cdf)  of  or  . 

4.1.4  APPLICATION  TO  NRL_LAND_CLUTTER  DATA 


FI6URE  4.4  CUMULATIVE  DISTRIBUTION  FUNCTION  OF  THE 
ESTIMATED  STANDARD  DEVIATION  OF  PHASE  ERROR  iNRL  DATA). 


Figure  4.4  shows  the  results  of  applying  the  above  procedure 
to  some  of  the  NRL  8-element  array  land  clutter  runs.  The 
transmit-receive  phased  array  operated  at  UHF  and  had 
spacing  between  elements.  The  transmitting  pattern  was 

Dol ph-Chebyshev  weighted  for  -£4  dB  sidelo! as.  The  system  was 
flown  repeatedly  over  the  south-eastern  portion  of  the  United 
States  at  about  £00  knots  at  an  altitude  of  15,000  to  £0, 000 
ft. 


Table  4.1  provides  important  information  relating  to  the 
family  of  curves  of  Figure  4.4.  The  curves  should  be  used  in 
conjunction  with  <4.  1).  For  example,  a  1  dB  loss  of  mainbearn 
gain  corresponds  to  «  5®0  mrad.  For  /Tl„  =  0.3,  95'/. 

of  the  candidate  bearnf orrners  result  in  a  loss  AG  <_  1  dB 
while  for  “’l  =  0.  15  a  ,11.  the  beam  formers  would  result  in  AG 

equal  to  or  better  than  0.6  dB. 


TABLE  4. 1  PERCENTAGE  OF  RANGE  BINS  IN  WHICH  THE  STANDARD  DEVIATION 

OF  THE  AMPLITUDE  DID  NOT  EXCEED  V 

A 


\  \ 

i 

.  15  1 

1 

1  ® 

1  OJ 

i  • 

1 

.  i£5  1 

.  30 

Np  j 

38  l 

1 

101  1 

1 

176  1 

1 

357 

nt  i 

1 

EH 

fii 

640  1 

1 

640 

Percentage  1 

1 

*/.5.  9  1 

1 

■xi  5.  a  i 

i 

*£7. 5  1 

1 

•X40 

N  =  Number  of  range  bins  that  passed  the  threshold  test 
N  =  Total  number  of  range  bins  searched 


Notice  that  the  results  of  Figure  4  confirm  the  expected 
rnonotonic  relationship  between  amplitude  variance  and  phase 
variance  upon  which  the  Minimum  (Amplitude)  Variance  algorithm 
is  based.  This  can  be  seen  in  the  change  of  slope  of  P(  q^)  as 
Tj,  varies.  Mathematically  speaking,  the  relationship  is 
embedded  in  the  Cauchy— Riernan  integrals  that  relate  the 
amplitude  and  phase  of  an  analytic  complex  function  (such  as 
the  electric  field  received  by  the  array). 

Notice  also  that  while  a  tighter  threshold  (smaller  ^  )  re¬ 

sults  in  a  set  of  beamformers  of  higher  quality,  the  resulting 
set  has  less  members.  The  tradeoff  between  these  two  factors 
(quality  vs.  number)  will  prove  to  be  important  while 
discussing  the  newly  developed  Robust  Minimum  Variance 
algorithms  later  in  this  report. 

0BE!=I£0II8N_IQ_iRIM_§QR_LAND_CLyTTER_DATA 

Although  the  results  of  the  NRL  data  are  highly  encouraging 
(for  example,  for  =  0. 15  about  6*  of  the  range  bins 

searched  contain  satisfactory  beamformers  (  4G(  0.6  dB)),  the 

results  seem  to  be  overopt imist ic  because  the  array  was  used 
both  on  transmission  and  reception.  It  is  simple  to  show  CEJ 
that  a  transmitting  aperture  of  size  D  will  result  in  a 
lobular  clutter  reradiation  patterns  of  a  nominal  lobewidth 
the  order  of  D  (in  the  vicinity  of  the  transmitter) .  This 
means  that  the  phase  of  the  received  field  is  expected  to 
slowly  vary  over  a  transmit-receive  aperture  and  our  chances 
of  finding  good  beamformers  are  good  to  begin  with.  A  more 
realistic  test  should  involve  a  bi static  situation  where  the 
transmitting  aperture  is  much  smaller  in  size  than  the 
receiving  aperture.  Although  we  could  not  find  clutter  data 
of  exactly  this  scenario,  we  realized,  however,  that  we  can 
use  SAR  data  in  a  way  that  fits  our  purposes. 


For  a  relatively  short  synthetic  array  of  a  length  much 
shorter  than  the  cross-range  dimension  of  the  illuminated 
patch,  the  transmitter  moves  over  a  small  distance  and 
illuminates  almost  the  same  patch  of  scatterers  at  each 
sampling  position.  The  field  sampled  by  the  synthetic  array 
will  be  almost  identical  to  that  received  by  a  bistatic  system 
having  a  receiving  array  the  same  size  of  the  synthetic 
aperture  and  a  transmitter  identical  to  the  SAR  transmitter. 

We  obtained  several  runs  of  land  SAR  data  from  ERIM  and  formed 
synthetic  arrays  of  a  maximum  length  of  Q  synthetic  elements 
out  of  a  total  of  10£4  available  sampling  positions.  Figure 
4.5  shows  the  results  of  a  typical  run  where  the  SAR  operated 
at  X-band  with  sampling  positions  taken  every  1.5m  (  a  50  “>•  >. 
The  transmitter  was  about  5  ft  in  the  cross-range  dimension 
and  had  about  -£0  dB  sidelobes. 

In  Figure  4.5,  we  set  =  0.  £  and  tried  different  synthetic 
array  lengths  <150^-  400^  ,  corresponding  to  3-B  synthetic 
elements).  As  expected  ,  for  a  given  the  chance  of 
finding  beamformers  (of  corresponding  quality)  decreases  as 
the  array  length  increases.  This  is  natural  because  it  takes 
a  small  (on  a  wavelength  scale)  dominant  scatterer  to  produce 
a  planar  phasefrorit  over  a  large  aperture.  Such  a  phase- 
front  will  be  certainly  good  for  phase-synchronising  a  smaller 
aperture  but  not  necessarily  so  for  a  larger  one.  Table  4. £ 
gives  the  percentage  of  range  bins  that  did  not  exceed  the 
threshold  relative  to  the  total  number  of  range  bins  examined. 


FIGURE  4.5  CUMULATIVE  DISTRIBUTION  FUNCTION  OF  THE 

ESTIMATED  STANDARD  DEVIATION  OF  OHASE  ERRORS  (SAR  DATA). 


There  were  1024  range  bins  per  run 


Whenever  we  needed  more  range 


TABLE  4.2  PERCENTAGE  OF  RANGE  BINS  IN  WHICH  THE  STANDARD  DEVIATION 

OF  THE  AMPLITUDE  DID  NOT  EXCEED  ^  =  0.2 


L  1 

1 

150  7 

i 

i 

200  7  1 

1 

250  7 

1 

1 

300  7. 

1 

1 

350  7- 

i 

i 

40071  1 

1 

Np  1 

1 

mm 

H 

277  1 

1 
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1 

1 

£5£ 

1 

1 

£38 

i 

i 

222  1 

1 

nt  1 

1, 024 

i 

i 

2,  048  1 

1 

4,  096 

1 

1 

5,  120 

1 

1 

9,  216 

i 

i 

14, 336 1 

1 

Percentage  1 

1 

£7% 

i 

i 

13.5V.  1 

1 

7.  4*/. 

1 

1 

4.  3% 

1 

1 

£  •  Et  */• 

i 

_!_ 

1 . 5%  1 

1 

bins  to  search  (to  get  a  reasonably  smooth  curve)  for  a  certain 
array  length  we  formed  additional  (nonoverlapping)  arrays  of  that 
particular  length.  We  had  no  trouble  doing  so  because  we  have  a 
total  of  1024  synthetic  array  elements  available  to  us. 

Figure  4.5  reveals  the  limitations  of  the  Minimum  Variance 
Algorithm.  The  algorithm  becomes  less  and  less  reliable  as  the 
array  length  increases.  For  example,  the  curve  for  L  =  400  7- 
(  si  12m,  about  the  si  zb  of  a  fighter  aircraft)  shows  that  if 

we  pick  a  single  range  bin  in  which  the  amplitude  variation  did 
not  exceed  'ylA  =0.2  (Table  4.2  shows  that  only  about  1.5V.  of  the 
available  range  bins  will  satisfy  this  threshold),  the  chance 

is  about  fifty-fifty  that  it  will  be  useless.  The  Minimum 

Variance  Algorithm  uses  a  single  range  bin  (the  one  that  has  the 
minimum  amplitude  variance  across  the  array)  as  the  reference 
signal.  Experiments  showed  that  this  range  bin  is  not  necessarily 
the  best,  and  for  larger  arrays  it  might  not  be  even  good  enough. 
As  we  will  show  later  (Under  Task  4.4),  we  developed  the  Robust 
Minimum  Variance  Algorithm  to  provide  high  quality,  reliable 
beamforming  based  on  a  £2ii.f=!££i.gn  of  inferior  beamf orrners.  For  L 
=  400  “X  (Figure  4.5),  we  were  able  to  reduce  05^  to  a  typical 

value  of  about  20  mrad.  The  algorithm  has  also  the  further 
advantage  of  allowing  for  subarray  processing  that  can  be  used  for 
sel f-coher ing  a  very  large  array. 

4.1.6  APPLICAIION_TO_SEA_CL.yiIER 

We  applied  the  Minimum  Variance  Algorithm  to  the  NRL  airborne 
Q-element  sea  clutter  data.  We  found  it  highly  unreliable. 
This  was  expected  because  of  the  obvious  low  chance  of  finding 
a  dominant  scatterer  in  a  sea  clutter  patch.  We  tried  it 
anyway  because  of  the  simplicity  of  the  algorithm  and  because 
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we  have  the  data  and  the  mechanism.  It  served,  however,  as  a 
check  for  the  principle  and  as  a  justification  for  going  to 
more  sophisticated  algorithms  for  self— cohering  on  sea 
clutter. 
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4.2  Task  2:  CLUTTER  DATA  SOURCES 

We  obtained  airborne  clutter  data  from  two  sources:  the  Naval 

Research  Laboratory  (NRL)  and  the  Environmental  Research  Institute 
of  Michigan  (ERIM).  We  are  quite  satisfied  with  the  amount, 
quality,  and'  diversity  of  the  data  that  we  acquired.  In  the 
following  we  give  a  brief  description  of  the  data  under  our 
disposal. 

4.2.1  NRL  8— Element  Array  Data 

The  Naval  Research  Laboratory  built  an  experimental  airborne 
phased  array  radar  system  that  was  flown  repeatedly  over  land  and 
water  and  collected  large  amounts  of  land  and  sea  clutter  data 
under  various  operating  conditions.  The  system  consisted  of  a 
rigid,  filled,  periodic,  3— element,  UHF  (  "X  s:  70  cm;  antenna  array 
that  was  used  for  both  radar  transmission  and  reception.  On 
transmission,  a  carefully  designed  Chybyshev  beam  was  launched 
(-24  dB  Sidelobes) ,  while  the  3  independent  receivers  were 
sampling  the  echoes  simultaneously  on  reception. 

In  each  receiver,  the  chirped  echoes  were  compressed  and  the 
coherently  detected  signal  from  each  range  cell  was  sampled  (with 
10-bit  precision)  and  stored  in  a  complex  in-phase  and  quadrature 
form.  The  range  cells  were  100  ft.  apart.  The  system  operated  at 
one  out  of  two  sub-KHz  pulse  repetition  rates  that  allowed  a  safe 
margin  of  oversampling  of  the  Doppler  spectrum  at  different  array 
pointing  angles  (aircraft  speed  was  held  to  about  200  knots  and 
the  altitude  varied  from  about  15,000  to  20,000  ft).  Array 
steering  was  mechanical.  The  number  of  pulses  (range  traces)  np 
and  the  number  of  range  cells  recorded  per  pulse  n*  was  chosen 
(according  to  memory  limitations)  such  that  np  :<  nR  a  1500. 

Land  clutter  data  runs  were  taken  from  flights  over  the  south 
eastern  portion  of  the  United  States.  We  have  runs  for  (np  ,  n„  ) 
=  (66, 192)  and  (np  ,  n*  )  =  (34,443).  For  most  of  the  runs  we 

obtained,  the  array  was  looking  along  the  ground  track  of  the 
aircraft.  We  requested  special  runs  for  which  the  array  looked  at 
about  45  degrees  relative  to  the  ground  track  to  allow  for  a 
reasonable  spread  in  the  Doppler  spectrum  to  test  the  Doppler 
Beam-Sharpening  algorithms  (Task  3). 

Sea  clutter  data  runs  were  taken  from  flights  over  the  Atlantic 
Ocean  off  the  south  eastern  shore  of  the  United  States.  The  data 
we  processed  corresponded  to  a  moderate  sea  with  the  aircraft 
flying  at  an  altitude  of  approximately  17,000  ft,  with  the  array 
broadside  pointing  at  the  direction  of  the  aircraft  nose.  For 
these  data  we  had  ( np  ,  n_  )  =  (3,5440). 
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4.£.£  ERIM  SAR  Data 

As  mentioned  earlier,  the  results  of  applying  the  minimum  variance 
algorithm  to  the  NRL  data  were  over  optimistic  because  the  array 
was  used  both  on  transmission  and  reception.  A  more  realistic 
test  should  involve  a  bistatic  situation  where  the  transmitting 
aperture  is  much  smaller  in  size  than  the  receiving  aperture. 
Although  we  could  not  find  clutter  data  exactly  fitting  this 
scenario,  we  realized  that  we  can  use  SAR  data  in  a  way  that  fits 
our  purposes. 

For  a  relatively  short  synthetic  array  of  a  length  much  shorter 
than  the  cross-range  dimension  of  the  illuminated  patch,  the 
transmitter  moves  over  a  small  distance  and  illuminates  almost  the 
same  patch  of  seatterers  at  each  sampling  position.  The  field 
sampled  by  the  synthetic  array  will  be  almost  indent  1 cal  to  that 
received  by  a  bistatic  system  having  a  receiving  array  the  same 
size  of  the  synthetic  aperture  and  a  transmitter  identical  to  the 
SAR  transmitter. 

We  obtained  several  runs  of  land  SAR  data  from  ERIM  and  formed 
nonover lapping  synthetic  arrays  of  a  maximum  length  of  S  synthetic 
elements  out  of  a  total  of  10£'4  available  sampling  positions.  The 
SAR  operated  at  x-band  with  sampling  positions  taken  every  1.5m 
(approximately  50  ).  The  transmitter  was  about  5  ft  in  the 
cross-range  dimension  and  had  about  -  £0  dB  sidelobes.  The  data 
were  intended  for  producing  SAR  images  (using  digital  processing) 
with  azimuth  resolution  of  about  1.44m  and  range  resolution  of 
about  1.5m.  The  data  were  compressed  in  range  and  each  range 
trace  consisted  of  10E4  range  cells. 
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4.3  Task  3:  DOPPLER  ALGORITHMS 
4.3.1  Doppler  Beam  Sharpening 

Figure  4.6  describes  an  airborne  radar  moving  with  speed  V  and 
illuminating  a  clutter  patch  with  beamwidth  £*■$  at  distance  R  and 
at  an  angle  9  from  the  ground  track.  The  receiving  array  of 
length  L  is  assumed  to  be  distributed  on  the  airframe.  For  a 
transmitter  of  size  D,  *  'X /l>  and  the  clutter  patch  width 

W— R  A0  «  MS/D  .  The  backscatter  from  the  clutter  patch  will  have 
a  lobular  pattern  of  a  nominal  lobe  width  (in  the  vicinity  of  the 
aircraft)  of  <  'X/kDR  «  D.  Typically  D<<L  and  such  a  received 
phasefront  is  useless  for  sel f-coher i ng  the  array.  By  virtue  of 
the  aircraft  motion,  however,  echoes  from  scatterers  within  the 
ground  patch  may  be  distinguished  from  one  another  by  their 
Doppler  shifts;  hence  narrowband  filtering  of  the  received 
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Effacttv«.  reduced  width 


FIGURE  4.6  NARROWBAND  FILTERING  IN  THE  RECEIVER  REDUCES 
EFFECTIVE  PATCH  WIDTH  (FROM  C31 ) 
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clutter  echoes  can  extract  the  reflections  from  scatterers  within 
a  subpatch  of  a  smaller  width.  The  smaller  the  subpatch  is  the 
wider  the  iobewidth  of  the  clutter  received  by  the  array.  To 
improve  the  chance  of  finding  a  planar  (or  spherical)  phase  front 
suitable  for  sel f-coher ing  the  array,  the  aircraft  must  be  covered 
by  a  relatively  small  portion  of  a  single  lobe  reradiated  from  the 
subpatch.  A  rule  of  thumb,  derived  in  C23,  requires  that  the 
subpatch  width  to  be  less  than  ^  R/2L  (corresponding  to  a 
Iobewidth  of  at  least  twice  the  size  of  the  aircraft).  The 
Doppler  shift  of  an  echo  from  a  scatterer  at  an  anglegfrom  the 
ground  track  is 


£  ft  2  vcasff/'X  (4.3) 

Using  (4.3),  it  is  shown  in  C23  that  a  narrowband  filter  of 
bandwidth  B  <  V  is  needed  to  pass  the  echoes  from  the  central 

“  L 

region  of  the  illuminated  ground  patch  of  Figure  4. 6  (  and  reject 
the  remainder.  This  requires  an  integration  time  Ta  B  )_  L.  The 

V 

NRL  data  that  we  Doppler-processed  were  collected  over  a  time 
interval  of  about  3L,  permitting  us  a  good  margin  for  testing 

V 

Doppler  Beam  Sharpening.  Figure  4.7  shows  the  results  of 
Dopplei — processing  an  NRL  data  run.  In  this  run  the  array 
broadside  made  45  degree  angle  with  the  ground  track  with  ( n* ,  n^ ) 
=  (66, 192) .  The  variable  parameter  here  is  Ns  ,  the  number  of 

pulses  (time  samples)  used  in  the  FFT.  For  each  curve  in  Figure 
4.7  an  Ns  -  point  Hamming  weighted  FFT  is  computed  for  every  array 
element  for  all  the  range  bins.  For  every  range  bin,  a  particular 
spectral  line,  correspond ing  to  a  certain  direction  in  the  clutter 


FIGURE  4.7  EFFECT  OF  DOPPLER  BERM  SHARPENING  ON  THE  CUMULATIVE 
DISTRIBUTION  FUNCTION  OF  THE  ESTIMATED  STANDARD 
DEVIATION  OF  PHASE  ERRORS  (NRL  LAND  CLUTTER  DATA) 
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patch,  is  examined  as  before.  That  is,  an  amplitude  variat 
threshold  is  set  on  the  normalized  RMS  of  the  magnitude  pattern 
that  spectral  line  across  the  array;  if  the  range  bin  passes 
threshold,  the  phase  profile  is  examined,  phase  unwrapping 
made-if  needed— fol lowed  by  linear  regression  and  an  estimate  of 


is  obtained, 
is  obtained 
d istr i but  ion 
range  bins 
total  number 


Repeating  the  procedure  for  all  range  bins,  a  set 
and  sorted  to  get  an  estimate  of  the  cumulat 
function  P(  ojp  )  .  Table  4.3  gives  the  percentage 
that  satisified  the  threshold  test  relative  to 
of  range  bins  examined. 
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TABLE  4.3  PERCENTAGE  OF  RANGE  BINS  THAT  SATISFIED  A  THRESHOLD  TEST 
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Notice  that  NT  is  always  a  multiple  of  =  19£.  In  order  to  get 
a  reasonably  smooth  curve  of  P<  ),  Np  is  increased  by  looking 
into  clutter  subpatches  of  different  directions  (correspond ing  to 
different  spectral  lines).  This  was  possible  for  Ns  >  £0.  For  Ns 
=  10,  in  addition  to  looking  at  two  different  directions,  we  had 
to  repeat  the  Doppler  processing  for  few  nonoverlapping,  well- 
separated  sequences  of  10  pulses.  For  Ns  =  1,  we  examined  several 
pulses  with  as  much  separation  in  time  as  possible. 


Notice  the  dramatic  improvement  in  the  statistics  of  ^ introduced 
by  Doppler  processing.  The  median  decreases  significantly  as  the 


number  of  pulses 
for  1  £  Nj  ^3 
slowing  is  due  t 
be  made  by  Dopp 


Most  of  the 


occurs 


pulses  Nj  increases.  Most  of  the  improvement  occurs 
N,  <  30.  It  slows  down  afterwards.  We  believe  that  the 
s  due  to  a  fundamental  limit  on  how  narrow  a  subpatch  can 
by  Doppler  processing.  It  is  implicit  in  (4.3)  that  the 


scatterers  are  assumed  to  be  completely  rigid  or  "frozen".  If  the 
scatterers  undergo  relative  motion  with  respect  to  each  other  due 
to  wind,  for  example,  as  in  the  case  of  vegetation  and  sea 
clutter,  the  Doppler  shift  due  to  a  particular  scatterer  will  not 
only  be  a  function  of  its  angular  position  but  also  a  function  of 
its  own  motion  relative  to  the  aircraft.  In  other  words,  clutter 
bandwidth  puts  a  limit  on  the  usefulness  of  Doppler 
beamsharpening.  While  it  is  evident  from  Figure  4.7  that  Doppler 
beamsharpening  is  very  useful  for  land  clutter,  it  is  shown  in 
C£3 ,  C3D  that  typical  clutter  bandwidths  of  sea  clutter  will 
render  it  useless. 


Notice  also  that  Figure  4.7  gives  a  unified  approach  to  the 
classification  of  Doppier  beamsharpening  algorithms.  The  Minimum 
Variance  algorithm  becomes  a  special  case  (Ns  =  1).  As  N 

increases  (Nj  >  1)  we  can  use  a  "Hybrid"  algorithm  where  we  use 

ail  the  pulses  available  to  us  in  a  Doppier  processing  that 
improves  our  chances  of  finding  dominant  scatterers  (by  reducing 
the  width  of  the  clutter  patch  and  hence  often  excluding  other 
scatterers  that  might  be  of  comparable  strength)  and  then  apply 
the  Minimum  Variance  algorithm  to  locate  (in  range)  such  a 
dominant  scatterer.  As  N$  becomes  large  enough  to  reach  the  region 
where  the  width  of  the  clutter  subpatch  cannot  be  made 
appreciably  smaller  (limited  by  the  clutter  bandwidth)  we  call  it 
a  "Doppler  Beamsharpening *'  algorithm. 

.  £  Minimum  Clutter  Bandwidth 

As  mentioned  earlier  in  the  previous  section,  the  above  techniques 
work  effectively  on  narrowband  clutter  (as  found  in  most  land 
clutter).  For  relatively  wideband  clutter,  as  in  the  case  of 
heavy  vegetation  and  sea  clutter,  we  have  two  techniques  (see 
Table  3.  1 )  : 

-1  Spatial  Correlation  (non-Doppler ) 

-2  Minimum  Clutter  Bandwidth  (Doppler) 

The  second  algorithm  C33  is  quite  complex  and  requires  a 
long  convengence  time.  It  uses  the  spectral  bandwidth  of  the  array 
output  as  a  "quality  measure"  to  determine  when  the  mainbearn  is 
best  formed.  It  was  developed  before  the  Spatial  Correlation 
algorithm  and  mainly  served  as  a  proof  for  the  existance  of 
algorithms  capable  of  self-cohering  on  wind  driven  land  and  sea 
clutter.  We  listed  this  algorithm  in  Table  3.1  for  theoretical 
interest  and  focused  our  attention  on  the  Spatial  Correlation 
algorithm  instead. 
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4.4  Task  4s  New  Algorithm  Development 


Introduct ion 


Me  developed  several 
study.  We  developed 
algorithm  to  provide 
inferior  beamformers  ( 
and  tested  the  Spatial 
to  sea  as  well  as  t 
processing  within  the 
cost  effective  large 
developed  a  techniqu 
above  new  algorithms  t 


new  algorithms  during  the  course  of  this 
and  tested  the  Robust  Minimum  Variance 
reliable  beamforming  on  a  collection  of 
targets  of  opportunity).  We  also  developed 
Correlation  algorithm  which  is  applicable 
o  land  clutter.  It  allows  also  for  local 
aperture  which  makes  it  attractive  for  a 
aperture  space  based  radar.  Finally,  we 
e  of  Subarray  Processing  for  applying  the 
o  very  large  airborne/spaceborne  radar. 


Robust  Minimum  Variance  Algorithm 


The  minimum  variance  algorithm  uses  a  single  range  bin  (the  one 
that  has  the  minimum  amplitude  variance  across  the  array)  as  the 
reference  signal.  Experiments  showed  that  this  range  bin  is  not 
necessarily  the  best  and  for  larger  arrays  it  might  not  be  even 
good  enough.  Also,  as  a  result  of  beamforming  on  a  single 
dominant  scatterer,  the  beamforming  angle  can  be  anywhere  in  the 
illuminated  sector.  Experiments  showed  that  once  an 
amplitude  variation  threshold  is  set,  a  corresponding  number  of 
range  bins  (on  the  average)  will  pass  the  test.  The  beamforming 
quality  of  these  range  bins  varies,  but  a  considerable  number  of 
these  will  have  quality  comparable  to  that  of  the  range  bin  of 
minimum  amplitude  variance.  The  idea  is  to  use  the  collection  of 
such  candidate  beamformers  in  a  phase-smoothing  operation  that 
results  in  a  superior  "equivalent"  phase  synchronizing  source. 


Figure  4. Q  illustrates  the  situation.  A  quasi— linear  receiving 
antenna  array  is  assumed.  The  array  is  distorted  both 
geometrically  and  electrically.  Let  us  assume  that  (M+l)  range 
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FIGURE  4.8  A  QUASI -LI NEAR  ANTENNA  ARRAY  IS  RECEIVING 
A  Range  Trace  That  Contains  M+l  Range  Bins  With  Dominant 
Scatterers. 
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range  bins  passed  a  low  amplitude  variation  t 
likely  that  each  of  them  contains  a  comma 
the  "degree  of  dominance"  might  vary  fr 
Assuming  far— field  conditions,  the  measured  p 
element  due  to  the  ith  range  bin  (that  passed 
have  the  general  form  (droppinn  anv  const a 
the  array) 


hold  " 1 1  .  It  1 
catterer  aithoug 
o  ne  t  o  a  n  o  trier 
by  the  nth  arra 
threshold)  wil 
chase  term  acres 


where 


<p.  = 

i,*i 


gives  the  angular  posi 
the  X-position  of  the 
determined  by  the  angu 
position  error  of  the 
deviation  from  planar 
unknown  phase  error  du 
assumed  to  be  uniform! 
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Notice  that  is  el 
will  be  almost  a  st 
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for  phase  unwrapping  by  adding  or  substract  ing  cl  H  -  After 
checking  (<£_<#»,•  )  for  i  =  and  phase  unwrapping  if 

necessary,  we  form  the  average 

V  a  .1  T* (<P  ~  =  ■£„  -  4r  *  2 * 

n  *  «l 

where  the  superscript  u  stands  for  "unwrapped"  and  the  bar  denotes 
averaging  with  respect  to  i.  We  get  a  smoothed  "equivalent" 
beamforrner  <p*  by  substract  ing  %  from  4>tn 


<&  =  4>n  -  % 


r»,n 

=  a,- 


S’  *«  >  +  *- 


J  n 


(4.  5) 


Comparing  (4.5)  and  (4.4),  we  recognize  (4.5)  as  an  equivalent 
phase  synchronizing  source  with  the  following  advantages: 
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var 
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The  approx irnately-equal  sign  accounts  for  the  fact  that  dominant 
scatterers,  found  in  nature,  are  not  necessarily  identically 
distributed  in  position  and  phase  deviations.  The  reduced 
variance  of  relates  directly  to  higher  beamforming  quality, 
while  the  reduced  variance  of  a,- means  a  reduced  uncertainty  in  the 
beamforming  angle,  i.e^,  the  formed  beam  will  have  the  extra 
advantage  of  being  centered  in  the  illuminated  sector. 


The  bj  term  in  (4.5)  should  be  treated  carefully.  Although 
geometrical  distortion  is  compensated  totally  along  the 

beam form ing  direction  in  the  simple  case  of  retrod  meet ive 

bearnf orrning  (as  applied  in  the  minimum  variance  algorithm),  a 
residual  geometrical  error  in  the  Y— dimension  will  be  present  at 
the  beamforming  angle  when  we  use  (4.5)  for  beamforming  (element 
position  errors  in  the  X-dimension  are  compensated  totally  at  the 
beamforming  angle).  The  reason  is  that  while  for  a  single 

beamforrner  a,=  ksin^.  ,  b  ,'  =  kcos  ft  ,  where  k  is  the  wavenumber  and 
9;  is  the  angular  position  of  the  beamforrner  relative  to  array 
broadside  (see  Apendix  B) ,  we  find  that  for  the  smoothed 

beamforrner,  the  beamforming  angle  9+  =  sin  {  a./k}  is  not 

necessarily  equal  to  cos  .  Fortunately,  we  will  show  that 

the  residual  geometrical  error  can  be  made  negligible  in  most 
practical  cases. 


PERFORMANCE 


At  the  beamf orming  angle,  the  response  to  a  test  target  will  be 
character i zed  by  a  residual  phase  error 


+  S>t(rO 

where 

Sfcn)  =  (b*  -  b,  )  %  , 

S<fi2(n)  =  K,n  , 

and  b#  =  kcos  .  The  standard  deviation  of  (4.6)  is 

°i+  «  < 

In  Appendix  B,  we  show  that 

t/a  _  l/2  _ _  x 

°s4>  *  15  £  { ( L  •  -  t  s.'ne.o1]  _  cos&i)  } 


(4.  6) 


(4.  7) 


(4.  Q) 


where  ,  is  the  standard  deviation  of  y 
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FIGURE  4. 10  GEOMETRY  USED  IN  THE  COMPUTER  SIMULATIONS 
USED  TO  ESTIMATE  er..  . 
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Although  we  were  able  to  derive  an  approximate  expression  for 
(4.8),  it  was  quite  complex  and  was  not  valid  for  small  values  of 
M.  We  prefered  to  estimate  through  computer  simulations 

where  we  assumed  that  the  angular  position  of  a  dominant  scatter 
is  uniformly  distributed  within  an  illuminated  sector  making 
an  angle  Q  with  array  broadside  and  having  a  width  &6  (Figure 
4.10).  Appendix  B  gives  the  details. 

Figure  4.11  plots  as  a  function  of  transmitter  beamwidth 

for  M=£0.  Gaussian  element  position  errors  in  the  y-di mens ion  are 
assumed  with  standard  deviation  of  5 y  .  The  transmitter  is 
assumed  to  make  a  10  degree  angle  with  array  broadside.  In 
practice,  a  typical  value  of  will  be  about  £  degrees  resulting 

in  a  negligible  geometrical  error  of  less  than  £  rnrad. 


FIGURE  4.11  GEOMETRICAL  PHASE  ERROR  AS  A  FUNCTION  OF 

TRANSMITTER  BEAMWIDTH  (COMPUTER  SIMULATION). 

Figure  4.  1£  plots  as  a  function  of  the  number  of  dominant 

scatterers  included  in  '  the  average,  M,  for  different  values  of 
transmitter  angle  &  .  A  beamwidth  &Q  of  £  degrees  is  assumed 

along  with  Gaussian  element  position  errors  of  =  5  'X  . 

Notice  that  the  standard  deviation  of  the  geometrical  phase  error 
vanishes  completely  for  M=l.  It  rises  sharply  for  1  <  M  £  3  and 

becomes  almost  independent  of  M  afterwards.  The  sharp  increase  of 
°it  as  &  approaches  S0  degrees  is  due  to  the  fact  that  then 
the  array  approaches  the  endfire  situation  where  errors  in  the 
y-direction  become  more  important  and  get  magnified.  Notice  that 
even  in  the  unlikely  scenario  of  &  —  80  degrees,  the  standard 

deviation  of  geometrical  phase  errors  ,  CJ^  ,  is  less  than  3  rnrad 
and  its  contribution  is  quite  negligible  relative  to  that  of  cr 
(see  (4.7)).  ^ 
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4100  'X  long  (see  Figure  4.5)  and  a  modified  version  of  the 
algorithm  was  used.  The  improved  version  of  the  algorithm 

selectively  includes  range  bins  that  already  passed  the  amplitude 
variation  threshold  according  to  whether  or  not  these  pass  a  phase 
threshold  set  on  the  RMS  of  the  regression  error.  Linear 
regression  is  done  after  phase  unwrapping  (Figure  4.3).  Provided 
that  the  dominant  scatterer  of  the  lowest  amplitude  variance  has 
sma 1 1  dev i at i ons  S«*<  n)  of  its  own,  a  larger  RMS  value  of  the 
regression  error  of  (  %  -  )  indicates  that  the  jth  range  bin 

should  be  excluded  from  the  average.  ft  low  RMS  value  of  the  same 
quantity  indicates  either  that  the  jth  range  bin  is  good  or  that 
£<n)  is  very  similar  to  ^  (n)  for  all  n  which  is  highly 

unlikely  for  an  array  of  large  number  of  elements.  For  an  array 
of  small  number  of  elements,  however,  one  can  exclude  such 

coincidences  by  not  relying  completely  on  one  range  bin  for 

testing  all  others.  More  than  one  range  bin  can  be  used  to  test 
the  rest,  and  a  majority  vote  or  a  linear  combination  rule  (or  the 
resulting  RMS  values)  can  be  used  to  selectively  include  or 
exclude  a  range  bin  (that  already  passed  the  amplitude  variation 
test)  from  the  average. 


Figure  4.13  shows  a  dramatic  reduction  in  oj.^  .  On  initial  value 
of  about  450  mrad  is  reduced  to  about  £0  mrad.  Assuming  a  typical 


An  initial  value 


residual  geometrical  error  of  5  mrad,  the  total  residual 
phase  error  <rs^  ss  (  =  &  mrad  which,  in  turr 
correpsonds  to  an  expected  loss  in  rnainbeam  gain  A  £  =  0.  00£  dB, 
which  is  quite  negligible  for  most  practical  cases. 


of  5  mrad,  the  total  residual 
) 1/1  =  £0.6  mrad  which,  in  turn 


In  conclusion,  the  robust  minimum  variance  algorithm  provides 
reliable,  high  quality  beamforming  using  a  collection  of  inferior 
beamformers.  In  addition,  the  formed  beam  is  centered  within  the 
illuminated  sector.  This,  by  itself,  is  an  important  advantage 
since  a  beam  formed  close  to  an  edge  of  the  illuminated  sector 
means  a  loss  of  almost  half  the  obtainable  field  of  view.  Also 
since  the  image  is  obtained  by  scanning  off  the  formed  beam, 
azimuth  measurements  will  be  relative  to  the  beamforming  angle. 
Since  the  formed  beam  is  now  centered  in  the  illuminated  sector, 
absolute  angle  measurements  will  be  much  more  accurate  (compared 
to  those  obtainable  by  applying  the  minimum  variance  algorithm). 
Finally,  and  as  we  will  show  later  in  this  section,  a  centered 
beam  allows  for  a  newly  developed  technique  of  subarray  processing 
that  can  be  used  to  self -cohere  a  large  distributed  array. 


As  a  final  remark,  it  is  worthwhile  to  notice  that  the  robust 
minimum  variance  algorithm  uses  the  range  bin  of  minimum  amplitude 
variance  in  an  interesting  way.  It  does  not  assume  that  it  is  the 
best.  It  only  assumes  that  it  is  good  enough  to  enable  phase 
unwrapping  of  other  range  bins  and  actually  have  a  fair  look  at 
their  phase  profiles.  The  robustness  comes  from  averaging  (or 
selectively  averaging)  these  phase  profiles. 
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4.4.2  SPATIAL  CORRELATION  ALGORITHM 


The  spatial  correlation  algorithm  enables  self-cohering  of 
distributed  arrays  in  the  absences  of  a  beamforming  point  source, 
using  the  spatial  correlation  properties  of  radar  clutter.  It  is 
shown  in  C43,  C53  that  if  the  spatial  autocorrelation  function  of 
the  field  (as  measured  by  adjacent  element  pairs)  is  ensured  to  be 
real  and  positive  in  the  neighborhood  of  the  origin,  both  periodic 
and  aperiodic  arrays  can  be  synchronized  forming  ret redirect ive 
beams  pointing  at  the  axis  of  synmetry  of  the  radar  transmitter, 
provided  that  the  interelement  spacing  does  not  exceed  some  limit 
(the  order  of  the  size  of  the  transmitting  aperture). 


Figure  4.14  shows  a  fairly  general,  quasi-linear  (geometrical ly 
distorted)  antenna  array.  The  array  can  be  both  thinned  and 
aperiodic.  The  elements  are  connected  to  a  signal  processor 
through  cables  or  any  other  sort  of  data  link.  The  quantities 
through  represent  a  random  phase  error  distribution  across  the 

aperture  due  to  uncertainties  in  cable  lengths,  impedance 
mismatches,  and  perhaps  a  slowly  varying  pattern  of  atmospheric 
turbulence.  In  addition,  we  will  lump  the  differential  phase  due 
to  the  angle  of  arrival  of  a  plane  wave  into  these  quantities. 
Thus,  if  we  can  estimate  the  S  ’s  we  should  be  able  to  form  a 
ret redirect ive  beam  by  substracting  them  out.  In  the  following, 
we  will  show  that  this  is  possible. 


/  i 

The  ith  element  delivers  a  measured  field  value  of  e,-  =  e,-  e 
where  e  is  the  error-free  value  (excluding  the  differential  phase 
due  to  angle  of  arrival  which  is  included  in  S,  )  produced  by  a 
large  number  of  non-coherent  scatterers.  The  ith  correlator 
provides  an  estimate  R»  of  the  quantity 


*  =  m  j  £  {« 


where  d,-  -  d  ,•  cos  ^  .  The  phase  of 

quantity  of  interest.  Assuming  that  we 
<  V*-  >,  we  can  approximately  write 


R'-  ,  called  V',- ,  is  the 

have  good  phase  estimates 


FIGURE  4. 14  ARRAY  MODEL  AND  SIGNAL  PROCESSOR 
CONFIGURATION  FOR  SELF-SYNCHRONIZATION. 
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yfif  s  arg  ^  •  Now  we  form  the  quantity 


P  '=  L  % 


(4. 10) 
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which  we  use  as  a  phase  correction  by  adding  it  in  the  ith  channel 
(for  all  i).  The  total  phase  shift  becomes 


l+<  *  +  t 


a  S. 


L  fi 

n»i  'n 


(4. 12) 


Now  the  phase  fiC.<)  of  the  spatial  correlation  function  plays  a 
major  role  in  synchronizing  the  array.  If  the  spatial  correlation 
function  is  real  and  positive  in  the  region  of  interest  (  /St  =  0 
for  all  i)  we  have  an  ideal  situation,  for  then  (4.1£)  reduces  to 


=  S'  *  constant  independent  of  t 


(4.  13) 


and  the  random  variable  has  been  replaced  by  a  constant.  The 
result  is  a  well-formed  retrod  irect  ive  beam  making  an  angle  Qa  with 
array  broadside. 


shown  in  C41,  C51  that  R(X)  is  the  Fourier  transform  of  the 


intensity  patterin  I(u)  associated  with  the  illuminated  clutter 
patch,  which  by  definition  is  real,  it  is  only  necessary  to  insure 
that  I (u)  be  symmetric  for  R(X)  to  be  real.  We  will  show  later 
that  this  is  a  practical  condition  to  achieve. 


PRACTICAL  REALIZATION 


The  technique  described  above  is  based  on  the  assumption  of  the 
availability  of  noncoherent  microwave  scatterers  and  on  our 
ability  to  accurately  measure  samples  of  the  spatial 
autocorrelation  function  of  the  reradiation  from  these  scatterers. 
A  large  number  of  nonrigid  scatterers  with  random  relative  motions 
can  provide  us  with  the  required  noncoherence.  Examples  are  sea 
waves,  moving  leaves  of  a  bush  or  a  field  of  scatterers  where  the 
size's  are  smaller  than  the  resolution  cell  of  the  observing 
instrument.  This  is  a  common  situation  in  radar.  We  can  estimate 


till 


R’  <X)  by  averaging  successive  correlation  readings  (cross 
products)  taken  at  time  intervals  larger  than  the  time  required  to 
decorrelate  the  geometry  of  the  individual  scatterers.  A  faster 
way  is  to  take  the  cross  products  from  echoes  corresponding  to 
successive  range  cells,  provided  that  the  geometrical  features  of 
the  scatterers  are  statistically  homogeneous  over  the  area  of 
interest.  Sea  clutter  is  a  good  example.  Even  a  rigid  but  rou  gh 
terrain  can  be  useful  if  echoes  from  successive  range  cells  are 
statistically  the  same,  i.e. ,  if  they  can  be  considered  as  sample 
functions  drawn  from  the  same  random  process. 

Figure  4. 15  shows  a  radar  transmitter  of  aperture  size  D.  It 
illuminates  a  far-field  patch  of  approximate  area  where 

is  the  transmitter  range  resolution  and  C>B  is  the  transmitter 
beamwidth.  For  an  observer  at  the  transmitter  site,  the  field 
produced  by  a  scatterer  at  point  1  is  equal  to  the  field  produced 
by  the  same  scatterer  moved  to  point  □  and  multiplied  by  the  phase 
delay  factor  e'1"1****"  .  Because  the  phase  relationship  between 

the  scatterers  at  points  1  and  0  is  completely  random,  the  phase 


FIGURE  4. 15  RADAR  ILLUMINATED  PATCH  AND  ITS  1-D  EQUIVALENT 

INTENSITY  DISTRIBUTION. 


difference  between  the  displaced  equivalent  of  the  scatterer  at  1 
and  the  actual  scatterer  at  0  will  be  random.  Thus  in  principle, 
all  the  scatterers  along  the  line  01  can  be  mapped  to  the  point  0 
and  the  whole  illuminated  patch  can  be  collapsed  into  an 
illuminated  arc  with  illumination  intensity  I(u).  If  the 
illuminated  patch  contains  a  large  number  of  statistically 
homogeneous  scatterers,  we  would  expect  the  intensity  I(u)  to  be 
proportional  to  the  square  of  the  radiation  pattern  of  the 
transmitter.  Thus  if  the  illumination  of  the  transmitter  is 
symmetrical-  so  too  will  be  I(u>,  and  R(X>  will  be  real.  For 
example,  the  radiation  pattern  of  a  uniformly  illuminated 
aperture  will  be  proportional  to  sincDu/'X  and  I<u)  will  be 
proportional  to  sine2  Du/^  .  The  autocorrelation  function  R(X), 
being  the  Fourier  transform  of  I(u>,  will  be  triangular  with  base 
2D.  In  general,  the  correlation  distance  will  be  the  order  of  D. 
Thus  the  largest  interelement  distance  in  the  receiving  array 
should  not  exceed  D  (with  some  margin).  Notice  that  the  plane 
wave  model  of  Figure  4. 14  will  be  a  good  approximation  for  most 
practical  situations  where,  typically,  the  transmitter  pattern  has 
a  sufficiently  narrow  beam  (few  degrees)  and  a  low  sidelobe  level. 

In  summary,  the  only  requirement  for  R(X>  to  be  ideal  is  that  it 
be  real,  which,  in  turn  requires  that  the  scatterers  be  uniformly 
distributed  within  the  clutter  patch  and  that  the  transmitter  beam 
be  symmetric  (which  is  typical).  A  real  taper  for  the 
transmitting  aperture  provides  the  desired  symmetry,  furthermore, 
from  the  symmetry  of  the  problem,  it  is  evident  that  the  beam  so 
formed  will  point  in  the  direction  of  the  axis  of  symmetry  of  the 
transrnitt  ing  antenna. 


£XPERIMENIS_WIIH_aiRBgRNE_SEa_CLUIIER_rtQD0R_D0in 

We  applied  our  algorithm  to  the  experimental  sea  clutter  data 
supplied  by  the  Naval  Research  Laboratory  (NRL).  The  8-element 
system  was  flown  over  the  Atlantic  Ocean  repeatedly  and  collected 
huge  amounts  of  sea  clutter  data  at  a  variety  of  operating 
conditions.  The  data  we  used  corresponded  to  a  moderate  sea  with 
the  aircraft  flying  at  an  altitude  of  approximately  17,000  ft. 
with  the  array  broadside  pointing  in  the  direction  of  the  aircraft 
ground  track. 

Because  the  array  was  rigid,  its  proper  performance  could  be 
predicted  and,  hence,  compared  to  the  performance  of  our  adaptive 
method.  Random,  unknown  phase  errors  uniformly  distributed  over 
C  ,  It  1  were  inserted  in  each  receiving  channel.  Figure  4.  16 
shows  the  results  of  applying  the  algorithm  to  one  such 
experiment.  This  excellent  agreement  is  typical  of  many  similar 
tests. 


When  clutter  patches  have  the  same  statist irs,  echoes  originating 
from  neighboring  range  ceils  can  be  treated  as  sample  functions 
drawn  from  the  same  random  process.  In  this  situation  samples  of 
the  spatial  correlation  function  R(X)  rnay  be  estimated  by 
averaging  a  large  number  of  correlation  readings  (cross  products) , 
each  taken  from  a  different  range  cell.  Thus,  echoes  due  to  a 
single  transmitted  pulse  were  used  to  phase  synchronize  the  entire 


FIGURE  4.  16  THE  RADIATION  PATTERN  OF  THE  6-ELEMENT  PERIODIC 

ARRAY  (FROM  L'SJ  )  . 

(  )  IDEAL.  (...)  DISTORTED  DUE  TO  PHASE  ERRORS. 

( - )  RECOVERED  BY  THE  SPATIAL  CORRELATION  ALGORITHM 

OPERATING  UPON  SEA  CLUTTER. 

array.  We  used  100  successive  range  cells  starting  at  6.4  miles 
and  ending  at  8.3  miles.  The  clutter  to  noise  ratio  was  about  30 
dB.  The  phase  of  the  estimated  autocorrelation  values  had  an 
estimated  standard  deviation  aj  of  0.097  rad  which  corresponds 
to  an  expected  mainbeam  gain  loss  of  0.089  dB  (from  (3.43),  L51 ) 

which  is  excellent  for  most  applications. 

A  small  beampointing  error  appears  in  some  experiments;  it  seems 
to  be  zero— mean  and  its  variance  is  believed  to  decrease  as  the 
number  of  antenna  elements  increases.  A  zero-mean  beampointing 
error  indicates  a  real  R(X).  This  agrees  with  what  we  expect 
since  the  transmitting  beam  is  symmetric  and  the  roughness  of  the 
ocean  is  fairly  homogenous  CGI. 


-?p_ 


cqncluding.remqrks 

In  conclusion,  we  think  that  it  is  worthwhile  to  emphasize  the 

following  important  features  of  the  spatial  correlation  algorithm. 

1.  Simple  and  fast:  it  requires  a  single  transmitted  pulse  for 
calculating  the  phase  error  corrections  necessary  for  real 
time  sel f— cohering  of  airborne  distributed  arrays. 

£.  Known  beamforming  angle  (if  the  clutter  is  homogeneous  as  in 
the  case  of  sea  clutter):  this  allows  for  the  application  of 
a  newly  developed  subarray  processing  technique  (described  in 
the  next  section)  that  can  be  used  to  phase  synchronise  a  very 
large  array. 


Suitable  for  sea  clutter  (as  well 


land  clutter). 


Allows  for  local  processing:  it  is  straight  forward  to  see 
from  (4.  10)  that  +  Vf  ,  which  means  that  the 
correct  ion  needed  for  the  next  element  is  equal  to  that  needed 
for  the  current  one  plus  the  angle  measured  by  the  correlator 
connecting  them.  This  means  that  each  element  need  only  to 
communicate  with  its  first  neighbors  (as  opposed  to  the  case 
of  the  minimum  variance  algorithm  where  the  amplitude  received 
by  each  element  should  be  sent  to  a  central  processor  to  test 
for  the  range  bin  of  minimum  amplitude  variance).  This 
property  makes  the  algorithm  particularly  attractive  for  a 
cost  effective  space  based  radar. 


v?,*r 


4.3  SUBARRAY  PROCESSING 


It  is  evident  from  previous  sections  that  the  larger  the  arr 
size  the  more  difficult  self-cohering  becomes.  Figure  i 
illustrates  the  point  for  the  minimum  variance  algorithm.  1 
robust  minimum  variance  algorithm,  however,  tolerates  lane 
arrays  as  long  as  the  phase  profile  of  the  range  bin  of  minin 
amplitude  variance  is  good  enough  to  permit  phase  unwrapping 
other  range  bins  that  passed  the  amplitude  variation  threshold  £ 
as  long  as  the  number  of  these  range  bins  is  large  enough 
effectively  smooth  the  phase  across  the  compensated  array.  It 
clear  that,  for  a  given  performance  specification  (and  oth 
design  parameters),  one  can  not  go  beyond  a  certain  maximum  art 


true  for  the  spatial  correlation  algorithm.  As  we 


showed  in  the  previous  section,  the  correction  needed  for  the  next 
array  element  is  equal  to  that  needed  for  the  current  one  plus  the 
angle  measured  by  the  correlator  correcting  them.  Thus  errors 
made  in  previous  correlation  measurement s  (in  angle)  will 
propagate  across  the  array  putting  an  upper  limit  on  array  size 
(for  a  given  performance  and  other  design  parameters).  The  effect 
of  error  propagation  is  analyzed  in  [51. 


It  is  logical  to  think  of  dividing  a  large  array  into  smaller 
subarrays  each  of  which  can  be  satisfactorily  self-cohered  using  a 
particular  algorithm  and  then  correcting  the  compensated  subarrays 
somehow  in  order  to  straighten  out  the  entire  array.  The  phase 


particular  algorithm  and  then  correcting  the  compensated  subarrays 
somehow  in  order  to  straighten  out  the  entire  array.  The  phase 
profile  across  a  compensated  subarray  will  be  approximately  a 
horizontal  line  with  an  unknown  height  that  d.ffers  from  subarray 
to  another  as  illustrated  in  Figure  4.17. 


'SUBARRAY 


^OVERLAPPING 

SUBARRAY 


FIGURE  4. 17  SUBARRAY  PROCESSING 
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The  subarrays  cannot  be  used  to  form  a  "superarray"  unless  phase 
Jumps  between  adjacent  subarrays  can  be  estimated  and  corrected 
for.  We  do  not  have  access  to  the  compensated  phase  profiles;  we 
only  have  access  to  the  corrections  that  our  sel f-coner ing 
algorithms  produced.  If  the  subarrays  are  made  to  overlap  over  a 
number  of  elements,  two  sets  of  corrections  will  be  calculated  for 
each  group  of  elements  shared  by  two  subarray  neighbors.  For  each 
such  element,  the  difference  between  the  two  corrections  provides 
an  estimate  of  the  phase  jump  between  the  two  compensated 
subarrays.  Naturally,  an  estimate  based  on  all  the  shared 
elements  will  be  of  higher  quality.  The  only  problem  with  making 
the  subarrays  overlap  is  that  we  are  making  them  longer  than 
necessary  this  way  (and  hence  making  the  effort  to  phase 
synchronise  each  of  them  properly  more  difficult).  A  better 
approach  is  to  let  the  subarrays  be  nonover lapping  and  form 
another  set  of  "overlapping"  subarrays  as  illustrated  in  Figure 
4.17.  The  phase  jump  between  two  adjacent  subarrays  can  be 
estimated  on  two  steps:  4.  a  (Figure  4.  17)  with  the  help  of 
an  overlapping  subarray. 

Notice  that  an  implied  condition  for  the  success  of  the  technique 
is  that  beams  formed  by  the  subarrays  should  be  pointing  basically 
in  the  same  direction.  This  condition  is  necessary  for  forming  a 
high  quality  "super beam"  based  on  the  entire  array.  This  almost 
excludes  the  minimum  variance  algorithm  as  a  possible  vehicle  for 
sel f-coher ing  the  individual  subarrays  because  then  the  subarrays 
are  not  necessarily  focused  on  the  same  dominant  scatterer  and  the 
formed  beams  will  be  looking  almost  anywhere  in  the  illuminated 
sector.  On  the  contrary,  the  robust  minimum  variance  and  the 
spatial  correlation  algorithms  fit  the  procedure  naturally  since 
all  the  formed  beams  will  then  be  centered  in  the  illuminated 
sector. 


CONCLUSIONS 


Airborne  distributed  arrays  can  be  reliably  self-cohered  on  either 
land  clutter  or  sea  clutter.  On  land  clutter,  several  algorithms 
exist.  Their  performances  can  be  improved  by  Doppler  processing 
and/or  robust  processing.  On  sea  clutter,  the  spatial  correlation 
algorithm  is  the  only  practical  technique  and  it  fits  the  problem 
nicely. 


Subarray  processing  is  necessary  for  very  large  arrays.  It 
employs  either  the  robust  minimum  variance  or  spatial  correlation 
algorithm.  It  makes  robust  processing  applicable  to  very  large 
arrays  and  limits  error  propagation  inherent  in  spatial 
correlation  processing. 


Finally,  local  processing,  achieved  by  the  spatial  correlation 
algorithm,  can  be  used  with  subarray  processing  for  a  cos1; 
effective  implementation  of  space  based  radar  where  the  bulk  of 
the  processing  is  done  Locally,  at  the  element  level,  using  only 
data  from  adjacent  elements.  This  obviates  the  need  for 
transferring  large  amounts  of  information  to  and  storing  it  in  a 
central  processing  node. 


! 
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6.  RECOMMENDATIONS  FOR  FUTURE  WORK 


L 


Work  should  be  continued  in  the  two  areas  of  spatial  correlation 
processing  and  robust  minimum  variance  algorithms  emphasizing 
subarray  processing  (as  a  means  of  self-cohering  very  large 
distributed  arrays)  and  space  based  3DI  applications.  For  both 
techniques,  consideration  should  be  given  to  (1)  an  initial  mode 
in  which  a  subarray  is  used  to  illuminate  a  larger  angular  sector 
and  the  full  array  is  used  as  multiple  beam  receiver  and  (£) 
iterative  sel f-coher ing  of  the  full  array  for  transmit /receive. 


6. 1  SPATIAL  CORRELATION  PROCESSING  STUDIES 


The  work  should  include,  but  not  be  limited  to,  the  following: 


6. i.l  Develop  a  two-dimensional  spat i a 1 -cor re 1 at l on  algorithm 
utilizing  subarray  processing.  This  involves  the  extension  of  tne 
one— d i mens ional  algorithm  to  a  two-dimensional  array  of  elements 
in  order  to  take  advantage  of  the  fact  that  the  required  phase 
corrections  can  be  estimated  more  accurately  when  correlations 
between  an  element  and  all  of  its  £-D  nei  gh  bors  are  utilized. 


6.1.2  Perform  analysis  and  computer  simulations  for  optimizing 
system  performance  in  terms  of  size  and  configurations  of  local 
clusters  of  elements  and  subarrays,  pulse  width,  number  of 
samples,  etc. 


6.1.3  Evaluate  the  mainbeam  gain  and  array  sidelobe  performance  in 
a  noise  and  jamming  environment. 


6.1.4  Establish  candidate  ^-D  self-cohering  systems  for  use  with 
active  sources  such  as  a  jammer. 


£  ROBUST  MINIMUM  VARIANCE  STUDIES 


6. £. 1  Conduct  optimization  studies  involving  subarray  processing 
of  large  arrays  with  the  robust  minimum  variance  algorithm  being 
used  for  self-cohering  the  individual  subarrays.  This  should 
include  investigations  about  the  optimum  array  size,  optimum 
amplitude  variation  threshold,  and  an  optimal  theory  for 
estimating  the  phase  jumps  between  subarrays. 


6. £. £  Develop  a  two-dimensional  robust  minimum  variance  algorithm. 
Phase  unwrapping  should  be  performed  over  tilted  planes  (instead 
of  the  1-D  straight  lines)  in  this  case.  Efficient  £-D  phase 
unwrapping  algorithms  should  be  investigated  (developed  if 
necessary ) . 


-37- 


£. 3. 3  Design  and/or  simulate  an  experiment  applying  the  robust 
minimum  variance  algorithm  in  a  simulated  SDI  threat  cloud 
environment  where  a  space-based  large  antenna  array  iteratively 
forms  a  t ransrni t /rece i ve  narrowbeam  on  a  dense  pack  of  flying 
space  object,  then  scans  the  beam  within  the  pack  in  order  bo 
image  the  threat  for  subsequent  identification  and  classification. 
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APPENDIX  A.  THE  RADIO  CAMERA  ALGORITHM 
A. 1  ADAPTIVE  BEAMFORMING 


exp(  j 


the 

3  (  i 

bin 


Table  A. 1  shows  all 
the  process:  A 
from  the  ith  range 
the  search  for  Re  . 
The  measured  phases  V? 

TfjC  +  HiT  .  The  first 
array  geometry  and  is 
of  array  distortion 
scattering.  The  second 


steps  in  the  procedure.  Step  1 
W*  )  is  the  complex  envelope  o 
received  by  the  nth  element. 
Step  3  is  the  adaptive  beamfor 
„  can  be  broken  into  the  sum  of 
term  contains  the  conventional 
all  that  would  be  expected  in  t 
,  medium  turbulence,  multi 


term  represents  the  errors  that 
TABLE  A. 1  STEPS  IN  RADIO  CAMERA  IMAGING  C 1 1 


i ni t i at es 
f  the  echo 
Step  £  is 
ming  step. 

two  terms 
target  and 
he  absence 
path,  and 
the  phase 


Step 

1  Measure  and  store  complex  en¬ 
velopes  of  echo  samples 


range  bin 


X  \ 


it. 


£  Find  Rtf  such  that  Aof,  «  A, 
all  n 


A  *■ 


element  number 

J 


Phase  rotate 
conjugate  in 


at  R,  by  phase 
relation  to  refei — 


ence  element,  exp  J  <  'VK,  -  Yi*  > 

4  Phase  rotate  at  all  range  ele¬ 
ment  s 

5  focus  at  each  range  R 

6  Phase  shift  linearly  with  angle 

7  Sum  at  each  range  element 


A  * 


i*. 


%o) 


■V" 


4v,« 


£  B,. 


S, 


-i  •<*..«* 


in 


S4(“> 
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nsi 

the  absence  of  error 
the  ith  range  element 


synchronisation  process  must  overcome.  In 
the  image  of  the  target  or  clutter  from 
would  be  obtained  from  the  integral  of  the  phase— weighted  signals 
received  across  the  array  from  the  ith  range  bin.  As  in  any 
phased  array,  the  phase  weighting  is  the  conjugate  of  the  kernel 
of  the  diffraction  integral.  For  ease  of  discussion  the 
diffraction  integral  can  be  approximated  by  the  Fresnel  integral, 
which  reduces  to  the  Fourier  integral  when  the  target  is  in  the 
far  field.  The  Fresnel  integral,  in  turn,  is  best  represented  by 
a  sum,  as  in  (A. 1),  because  the  array  is  discrete.  The  sum 

*  £  i  //„  -jkf*na  -Zn/2RA 

Sj‘ tu)  —  A  Vi  ^  ^ 


(A.  1  ) 


is  the  image  s(-  of  the  scattering  sources  in  the  ith  range  bin, 
provided  that  y"  -  0.  We  have  u  =  sin©  and  9  —  scan  angle 
from  the  normal  to  the  array.  The  discrete  variables  x„  and  Rt- 
are,  respect ively,  the  x-coordinate  in  the  array  (transverse  to 
the  array  normal)  of  the  nth  element  and  the  distance  to  the  ith 
range  bin. 

* 

The  function  of  the  adaptive  processor  is  to  compensate  for  VLto 
allow  the  operation  described  by  (ft. 1)  to  be  completed.  In  the 
earlier  steps  of  the  process  the  data  were  searched  to  find  the 
range  bin  in  which  the  echo  amplitude  across  the  array  had  the 
smallest  variation.  That  range  bin  became  the  reference  range 
and  its  signal  became  the  synchronizing  signal.  Its  complex 
envelope  is  ft0„exp(  j  V'o „  ),  where  the  amplitudes  are  all  nearly 
the  same.  The  phases  of  the  echoes,  however',  are  random  because 
of  the  perturbed  geometry  of  the  array,  the  spatial  variations  of 
the  refractive  index,  the  impedance  variations  from  element  to 
element,  or  any  combination  of  these  factors.  If  the  reference 
source  were  ideal,  the  phase  differences  would  be  due  entirely  to 
them.  In  addition,  the  phases  are  perturbed  by  multipath  and 
scattering.  But  because  it  is  generally  impossible  for  the 
system  designer  to  obtain  a  priori  information  about  these 
conditions  the  signal  processor  must  necessarily  ignore  them. 

The  third  step  in  the  process  is  to  compensate  for  the  phase 
variations,  which  are  assumed  to  be  due  to  the  first  set  of 
factors  (i.e. ,  geometry,  refractive  index,  and  impedance; . 
Correction  is  made  by  phase-rotat ing  the  complex  envelopes  of  the 
signals  from  Ra  received  by  the  different  array  elements  to 
cophase  them.  Let  some  element  in  the  array  be  designated  the 
reference  element  and  the  phase  of  its  signal  from  the  reference 
range  be  designated  “Vi,  .  The  proper  phase  shift  for  the  nth 
element  is  the  negative  of  the  phase  difference  ■■=  Yon  -  'Vi* 

The  complex  signal  envelope  at  the  nth  element  becomes  ft  0„  exp 
(  j  Vi„  )expC-  i  (  -  Vit  )  1  Si  ft  exp(j  >.  This  correction,  or 
phase  conjugation,  is  exactly  what  a  phased  array  or  lens  would 
do  when  focusing  on  the  reference  source.  The  output  signal  from 
the  array  when  it  is  focused  is  the  sum  of  these  phase-correct ed 
echoes  from  Ra  . 


FIGURE  ft. 1  COMPONENTS  OF  THE  PHftSE  SHIFT  IN  ftDfiPTIVE  BEftMFORM ING. 


'the  differential  distance  •'•hat  it  represents,  is 

the  sum  of  two  quantities.  The  first  is  due  to  the  geometric 
distortion  and  includes  the  steering  direction  to  the  source.  It 
is  the  distance  of  each  element  from  a  datum  line  or  plane  in  the 
array  that  contains  the  reference  element  and  that  is  oriented 
normal  to  the  direction  to  the  synchronizing  source,  which  is  the 
desired  correction.  When  the  source  is  in  the  near  field  of  the 
array  also  has  a  second  component.  it  is  the  added 

differential  distance  »,f/£R  between  the  spherical  phasefront  and 
the  datum  line.  These  components  of  are  illustrated  in 

Fi gure  ft. 1 . 

ft. £  SCANNING 

The  fourth  step  (performed  simultaneously  with  Step  3)  is  to 
phase-rotate  all  the  signal  samples  from  each  antenna  element  ay 
The  samples  of  the  complex  envelope  from  the  ith 
range  bin  now  become  A  exp  L  j(  -  Y£„  +  )  1 .  The  sum  of  a 

set  of  complex  samples  represents  the  output  of  a  misfocusea 
phased  array  or  lens  because  ^  ,  except  at  the  range  of 

the  reference  reflector. 

The  fifth  step  is  to  focus  the  array  at  all  ranged 
simultaneously.  This  task  is  accomplished  for  an  arbitrary  range 
R,- by  refocusing  the  array  from  the  reference  range  R#to  R,-.  The 
phase  correction  is  approx imately  quadratic,  given  by  (A.  1). 
Assuming  that  the  earlier  steps  were  performed  properly,  the 
self-cohering  process  forced  the  quadratic  component  of  the  phase 
shift  of  the  signal  in  the  nth  channel  to  become  -kx*/£Rt-  + 
k  xV£R0 .  To  focus  the  array  to  range  Rj  this  term  must  be  set  to 
zero,  which  requires  a  further  phase  addition  of  (kx£  /£)  (1/R{- 
1/R,  ).  This  step  requires  a  knowledge  of  the  range  R  of  the 
reference  reflector.  Fortunately  the  value  of  R0  can  be  measured 
(as  in  conventional  radar)  by  the  round-trip  travel  time  of  the 
pulse  to  the  phase  synchronising  source  and  it  is  read  directly 
into  the  signal  processor  from  the  radar  receiver.  The  accuracy 
of  measurement  is  determined  by  the  range  resolution  of  the 
system,  which  is  the  order  of  the  reciprocal  of  the  signal 
bandwidth  .Cin  distance  units]  or  the  near-field  range  bearnwidth, 
whichever  is  smaller. 


-A4- 


Step  6  imparts  a  linear  phase  rotation  to  the  range— f ocusea 
complex  envelope  (designated  B,-*  in  Table  ft.  1)  for  each  scan 
angle  u.  The  phase  shift  is  -kxn  u. 

The  last  step  forms  the  sum  of  the  linearly  phase-weighted, 
range-focused  samples  to  obtain  the  image  si (u>  for  the  ith  range 
bin : 


s =  Z'  6,,  t 


-i 


ft.  3  ELEMENT  POSITION  TOLERANCE 


(ft.  £) 


Steps  £,  3,  and  4  require  no  knowledge  whatsoever'  of  element 
position.  The  phase  synchronization  process  is  purely 
ret rod i rect i ve.  Steps  5  (refocusing  in  range)  and  S  (scanning  in 
angle)  do  require  coordinate  information.  The  tolerances  on 
element  position  error  have  been  worked  out  C1J,  LEI.  The  most 
stringent  tolerance  is  invoked  by  Step  6.  Briefly,  the  theory 
consists  of  the  following  points: 


First,  the  loss  in  mainlobe  gain,  in  dB,  due  to  all  the  random 
phase  errors  across  the  array  is  shown  to  be  ~  4.  3  <r$  ,  where 
is  the  variance  of  the  phase  errors  in  square  radians.’ 


*■; 


Second,  the  phase  variance  due  to  random  position  errors  is 
approximately  «  kx »£ *  ,  where  OJ*  is  the  position  error 
variance  in  the  array  in  the  direction  perpendicular  to  the 
beamforming  direction  and  Q  (rad)  is  the  scan  angle  measured 
from  the  direction  of  beamforming.  Combining  these  expressions 
and  allowing  1-dB  loss  in  array  gain  as  an  accepatable  gain-loss 


tolerance,  the  rrns  element  tolerance  becomes 


*  (  V/4 TT  )  e\ 


at  the  radar,  &, 


The  factor  'X  /47r  is  the  conventional  tolerance  in  phased  array, 

mirror,  or  lens  design;  Q is  the  maximum  scan  angle  from  the 
direction  of  the  reference  source  and  is  half  the  field  of  view. 
Because  an  individual  radar  target  always  subtends  a  small  angle 

(<  1  and  the  allowed  element  position  error  is 
for  example,  if  the  system  were  designed  to 
large  (in  angle)  as  the  moon,  $mtki  would  be 
rad.  In  this  case  the  position  accuracy 
by  two  orders  of  magnitude  to  about  .  It 

extraordinary  liberty  in  position  tolerance, 
phase  synchroni zat ion,  that  permits  radio 
a  nonrigid  or  ill— surveyed  array. 


exceedingly  large; 
image  a  target  as 
approximately  10"2 
tolerance  increases 
is  precisely  the 
following  adaptive 
camera  imaging  with 
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APPENDIX  B.  STATISTICS  OF  GEOMETRICAL  PHASE  ERRORS 


FIGURE  B. 1  A  TWO-DIMENSIONAL  GEOMETRY  OF  A  QUASI— LINEAR 
ARRAY  EXCITED  BY  A  POINT  SOURCE. 

Figure  B. 1  illustrates  a  quasi-1 inear  array  laid  along  the  x-axis 
and  excited  by  a  point  source  making  an  angle  6  with  array 
broadside.  From  the  geometry,  the  distance  &  between  the  source 
at  (x,  ,  y(  )  and  the  nth  array  element  at  (x„  ,  y„  )  is  given  by 

f„  =  t  +  (< j, 

=  [  CR  y.Vifl  .  Xn)x  +  (Reas*  .  *„)*  yn 

=  R  t 1  ♦  £  -  (  xn  s.'"0  -*■  co\9)]'>L  (B.  i  ) 

Expanding  (B. 1)  and  retaining  the  first  order  terms  in  the 
binominal  expansion  we  get  the  Fresnel  (near  field)  approximation 

/»  ~  R  +  'OS*)  (B.L) 

If  we  further  assume  that  <r  /R>  <<1,  we  get  the  Fraunhofer  (far 

field)  approximation 

f„  ~  ?  -  (Xn  S«Vi0  +  >jn  cos  9)  (B.  3) 

and  the  difference  in  phase  retardat ion  (relative  to  the  origin) 
sensed  by  the  nth  element  is 

—  If  ^  f„  s  If  in  -V  If  COS  &  'if)  (B.4) 


If  we  have  a  number  of  non-ideal  isolated  sources  (separated  in 
range),  the  phase  measured  at  the  output  of  the  nth  receiving 
channel  of  the  distorted  array  due  to  the  ith  source  is 


s  *«•  +  k  coift  ~jn  +  S 

»#W  2tt 

-  +  b>  %  +  +€» 

where 


(B.  5 > 


m 


t> 


ksin  0}-  , 
kcos  , 

=  deviation  from  planar  phase  front, 
unknown  phase  error  pattern  due  to  array  electrical 
distortion;  assumed  to  be  uniformly  distributed  be¬ 
tween  —  TT  ,  -t-  IT  , 

integer  to  account  for  the  modulo  £  nature  of  phase 
measurements. 


As  discussed  in  section  (4.4)  if  we  have  (M+l)  sources,  the  one 
that  has  the  minimum  amplitude  variance  is  used  to  phase  unwra 
the  rest  so  that  their  phases  can  be  averaged  to  produce  th 
equivalent  smoothed  source 


4$ *  s  £>'  +  k '  *in  ♦  +  €, »  -t-nr  i„  (4.5) 

In  the  radio  camera  algorithm  (Appendix  A),  the  beamformer  phas 
is  substracted  from  the  phases  of  all  range  bins  and  a  scannin 
term  is  applied  to  image  each  range  bin  in  azimuth.  To  evaluat 
the  quality  of  (4.5)  as  a  beamformer,  we  first  substract  it  from 
the  array  response  to  an  ideal  test  target  (a  point  source) 

4>n  s  Qr  X*  +  +  2n  /Tn  (B.  6) 

The  difference  (modulo  £  V  )  is 


%rs  ~  +*  -  (  *T  ~  ■+  C  bT  -*)*„  -  <B.  7) 

We  then  scan  the  test  target  in  angle.  The  peak  of  the  magnitude 
pattern  occurs  (within  a  very  good  approximation)  at 


a.  m  m  cn  tu 


and  the  residual  phase  error  at  the  beamforming  angle  becomes 

r*  H 

ss  S^(n)  +  (B.Q) 

where 

=  a,-£)%  » 

***<«>  =  - 

and 

b^  =  It  cos  {  Sin  ( 0;  /  It)  ^  . 

Assuming  that  the  element  position  error  in  the  y-d imension,  4n  i 
is  a  zero  mean  random  variable  independent  of  bt-  and  i  it 

is  straighforward  to  show  that  the  covariance  between  Sf((n)  and 
S^(n)  is  zero.  Thus,  the  standard  deviation  of  (B.8)  is 


<%  =  [  <r^  +  <rs\  ]V1  (B.s) 

In  the  following,  we  simulate  S4,  (n)  (and  hence  estimate  > 

while  we  present  typical  values  of  in  section  4.4  (based  on 

ERIM  SAR  data).  '  1 

The  variance  of  (n)  is  calculated  as 

iT'g.  =  v*r  [  t  b*  -  b.  ) 

■=  c  KvW'IeK]- 

(B.  10) 


where  we  took  advantage  of  the  assumed  independence  of  y0  and 
both  b*,  b,-  .  Furthermore,  since  y^  is  assumed  to  be  zero-mean, 
(B.  lid)  is  simplified  to 


That  is 


=  ®  |  ( b#  -  b.' )  } 

=  ir’  rvl  £  {( cos  C  s;«'(  s.«  -  cos©,)1  } 

%t  =  «  %  Edit'd)} 


(B.  1  i) 


(B.  1£) 


where 
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% 
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-f  ($)  =  C  I  -  C  Si'n  5,0  ^  -  C0S  9i 


(B. 13) 


UJe  estimated  E  <f  (  §  )  >  by  simulating  f(  5)  on  the  computer  where 
we  assumed  that  the  angular  position  of  a  dominant  scatterer  9;  is 
uniformly  distributed  within  an  illuminated  sector  making  an  angle 
9  with  array  broadside  and  having  a  width  &&  (Figure  4.10). 
Equation  (B. 12)  is  then  used  to  estimate  .  In  all  the 
simulation  experiments  we  performed  the  sample  'size  was  always 


made  large  enough  to  ensure  smooth 


curves. 


■Mill 


